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Abstract 

We consider structure formation and cosmic microwave background (CMB) aniso- 
tropics in a closed universe, both with and without a cosmological constant. The CMB 
angular power spectrum and the matter transfer function are presented, along with a 
discussion of their relative normalization. This represents the first full numerical evolu- 
tion of density perturbations and anisotropics in a spherical geometry. We extend the 
likelihood function vs. from the COBE 2-year data to > 1. For large the presence 
of a very steep rise in the spectrum towards low £ allows us to put an upper limit of 
< 1.5 (95%CL) for primordial spectra with n < 1. This compares favorably with 
existing limits on 0. We show that there are a range of closed models which are con- 
sistent with observational constraints while being even older than the currently popular 
flat models with a cosmological constant. Future constraints from degree scale CMB 
data may soon probe this region of parameter space. A derivation of the perturbed 
Einstein, fluid and Boltzmann equations for open and closed geometries is presented in 
an appendix. 

Subject headings: cosmology:theory — cosmic microwave background — large-scale 
structure of universe 



1 Introduction 

Much attention has recently been paid to structure 
formation and cosmic microwave background (CMB) 
anisotropics in open universes (e.g. Sugiyama k Silk 
1994, Kamionkowski k Spergel 1994, Kamionkowski, 
Spergel k Sugiyama 1994, Kamionkowski et al. 1994, 
Hu, Bunn k Sugiyama 1995, White k Bunn 1995, 
Gorski et al. 1995, Cayon et al. 1995). The idea 
that < 1 has been seen as perhaps more promis- 
ing because of the possibility of fluctuation genera- 
tion in open inflationary models (e.g. Lyth k Stew- 
art 1990, Ratra k Peebles 1994, Bucher, Goldhaber 
k Turok 1994, Lyth k Woszczyna 1995, Yamamoto, 
Sasaki k Tanaka 1995, Liddle et al. 1995). In this pa- 
per we would like to focus instead on the alternative 
case of a universe with a spherical spatial geometry, 
i.e. a closed universe. We will apply some of the re- 
sults and techniques developed recently in the study 
of open and flat spatial geometries to the closed case. 

Historically, there has been much interest in closed 
universes (Wheeler 1968, Hawking 1984), since the 
compact spatial surfaces can make quantum field the- 
ory more clearly consistent than in hyperbolic spaces 
and renders tractable the integrals of quantum grav- 
ity. The idea of the need for closure goes back to Ein- 
stein (1934), who regarded it as necessary to solve the 
"problem of inertia" , and was postulated as a bound- 
ary condition by early workers in the field (see discus- 
sion in Misner, Thorne k Wheeler 1970, p. 543, 704). 
Recently Linde (1995) has shown that it is possible 
to produce inflationary models which result in closed 
universes, in which the universe is "created from noth- 
ing" (Tryon 1973, Zel'dovich 1981, Zel'dovich k Gr- 
ishchuk 1984). This is possible because a closed 
universe has zero total energy, just as it has zero 
total momentum and total charge (Landau k Lif- 
shitz 1975). An extensive list of references to the 
early literature on closed universes can be found in 
Bjornsson k Gudmundsson (1995). 

The cosmological constant is often reconsidered in 
times of apparent astrophysical crisis (see Carroll, 
Press k Turner 1992), and today's Hubble constant 
versus age discrepancy is no exception. Recent in- 
dications that Hq ~ 80 kms _1 Mpc _1 partially moti- 
vated us to consider models with cosmological con- 
stant A > 0. As we discuss in §4, it is possible to 
gain a few Gyrs in such models, although exactly how 
old the universe can be depends largely on gravita- 
tional lens constraints (and also on constraints from 
large-angle and degree-scale CMB anisotropics). But 
certainly A is not the universal panacea for the age 
problem. 



While we allow for a non-zero A, we restrict our 
attention to those models which start at arbitrar- 
ily small scale (the Big Bang). Models without 
a Big Bang (i.e. having excessive amounts of A) 
can be excluded on quite general grounds (Ehlers k 
Rindler 1989). In order to specify an initial spec- 
trum of fluctuations, we shall implicitly assume an 
early period of inflation which gave rise to scale in- 
variant potential fluctuations (though we do not cal- 
culate the spectrum arising from any particular in- 
flationary model). Note that we only concern our- 
selves with scalar perturbations here, neglecting the 
possibility of tensor modes. We further restrict our 
attention to the well-motivated case of adiabatic den- 
sity fluctuations, although many of our general results 
would apply to isocurvature fluctuations as well, and 
it would certainly be possible to perform calculations 
for any specific isocurvature spectrum. We also con- 
sider what would happen if the last scattering surface 
was at the geometrical antipode (a putative "solu- 
tion" of the smoothness problem), although in reality 
this situation is inconsistent with current data. 

The outline of the paper is as follows. The next 
section reviews the closed universe geometry and es- 
tablishes our notation. §§3 and 4 present some mate- 
rial which will be necessary to understand our princi- 
ple results. The rest of the paper presents the results 
of our numerical calculations. §5 discusses the CMB 
fluctuations in a closed universe. Here we perform 
the first detailed calculations of the anisotropy an- 
gular power spectrum in closed models, by accurate 
integration of the coupled Einstein, fluid and Boltz- 
mann equations (for more detials see Hu et al. 1995, 
hereafter HSSW, and references therein). An exten- 
sive appendix gives a derivation of these equations for 
arbitrary geometry. §6 reviews the relation to large- 
scale structure, with explicit calculation of the matter 
transfer function and accurate calculation of the nor- 
malization relative to the CMB, including the effects 
of growth rate and non-trivial gravitational potential 
evolution. We show that the continuation to closed 
geometries is benign. In §7 we use the COBE 2-year 
data to place an upper limit on f2tot> and §8 reviews 
an extreme model suggested by Harrison (1993). §9 
reviews other limits on closed universes, and §10 con- 
tains our conclusions. 

2 Spherical Geometry 

We will use the following notation for contribu- 
tions to the critical density /9 cr ; t = 3Hq/(8ttG): fl ra , 
£l r , Q\, fix, arising from the (pressure free) matter, 
the radiation (photons plus massless neutrinos), the 
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cosmological constant and the curvature, respectively. 
The Friedmann equation can be used to relate these 
fi's: 

1 = n m o + fii-0 + ^AO + &K0- (1) 

We will generally use these quantities to parame- 
terize the present state of the universe, and hence 
the subscript '0' will be implicit throughout the pa- 
per. We will also sometimes refer to the quantity 
fltot = &m + fi r + f^A- The curvature contribution 
can then be written as fix = 1 — fitot- Note that for 
a closed universe fix < 1. 

The closed universe has metric 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



Kr 2 



(d0 2 



>d<f> 2 



(2) 



1) is the "curvature", and 



where K = -ffo(fi tc 
we use units where c = 1. Ho is the Hubble con- 
stant today, which we will usually denote by h = 
iZo/lOOkms^Mpc" 1 . A closed universe has curva- 
ture K > 1, i.e. a real radius of curvature. The 
open universe can be obtained by analytic continu- 
ation onto the imaginary axis. 

The metric can also be written in terms of the 
development angle x an d the conformal time r\ = 
J dt/a(t), as 

ds 2 = a 2 (r])[-dr] 2 + K' 1 {d X 2 + 

sin 2 * (d0 2 + sin 2 6><^ 2 )}] . (3) 

Intervals along geodesies are then related by d\ = 
\[Kdr\. Note that in terms of these coordinates (3) is 
just the metric of a 3-sphere, with us at the "North" 
pole. It is sometimes useful to imagine this geometry 
as a sphere embedded in a 4-dimensional Euclidean 
space (w, x, y, z) by (Weinberg 1972) 



w 
x 

y 



cos\ 

sin x cos 9 
sin x sin 9 cos c 
sin x sin 9 sin (, 



(4) 



In terms of this picture all null geodesies (light 
rays) form great circles on this sphere, with a max- 
imum separation when x = tt/2. A particularly im- 
portant point in the spherical space is the antipode, 
X = 7r, where the extrapolation of any two geodesies 
arriving at our position will cross. Thus for light rays 
with fixed angular separation, the coordinate distance 
subtended increases from zero (at * = 0) to a max- 
imum (at x = 7r /2) and then returns to zero again 



(at x = 7r )- Equivalently, an object placed at the an- 
tipode will fill the whole sky. We will see that this 
structure has important consequences for the angular 
power spectrum of CMB anisotropics. 

3 Inflation, Evolution 8z Flatness 

In a closed universe, with A = for the moment, 

2 



1 + cos fj 



(5) 



where f\ = \JKr\ is the development angle, which runs 
from to 27r. This means that if you live in such 
a universe, you can tell what stage of evolution you 
are at by measuring fi ra alone. The universe evolves 
from fi = 1 at the Big Bang, to fi — > oo at maximum 
expansion, when /9 cr ; t oc H 2 — > 0, and then back to 
fi = 1 at the Big Crunch. 

One consequence of this is that the universe spends 
much of its time near fi = 1. For these A = mod- 
els then, an argument can be made (A. Linde, private 
communication) which is a twist on the old Dicke flat- 
ness argument (Dicke 1970, Dicke & Peebles 1979). In 
a closed model with A = 0, the universe has 1 < fi < 2 
for half of its total conformal time (this has to be re- 
duced by another factor of 2 if you are only interested 
in the expanding phase). [In terms of cosmic time 
t rather than conformal time r/, this is about 18%, 
or half of that if the contacting phase is neglected.] 
Hence it is relatively likely that we would find fi not 
far from 1 in a closed model. It could then be argued 
that if we observe fi ~ 1, there is no flatness problem 
providing that fi > 1, and hence we probably live in a 
closed universe. Of course these are just a posteriori 
arguments based on one sample, and therefore may be 
entirely meaningless (see e.g. Evrard & Coles 1995). 

On the other hand, as is well known, the flatness 
problem is exacerbated when A ^ (e.g. Martel & 
Wasserman 1990), since the value of fi ra deviates very 
rapidly from unity near the present time, when A 
dominates the expansion. Thus there is a natural- 
ness problem (why is fi ra ~ ~ 1?) even in flat 
(i.e. fi ra + fiA = 1) models with A ^ 0. We shall con- 
sider models with fiA > despite these philosophical 
problems (Weinberg 1989), taking the view that even 
unlikely regions of possible parameter space deserve 
some exploration. 

We shall implicitly assume throughout this pa- 
per that the initial spectrum of density perturba- 
tions is scale invariant (in the gravitational poten- 
tial, see later) and sufficiently small to not overpro- 
duce fluctuations in the CMB. Such assumptions can 
be motivated by inflationary models which result in 
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fltot > 1- Self-consistently we will neglect any con- 
tribution to the CMB anisotropy from gravitational 
waves, which should be sub-dominant for scale invari- 
ant spectra. Since "standard" inflation generically 
gives fltot = 1) but allowing for a A component (Pee- 
bles 1984, Turner, Steigman k Krauss 1984), we will 
assume that we can have A in some fl t ot > 1 models 
also. 

Recently Linde (1995) has shown that it is possible 
to produce inflationary models which result in closed 
universes, in which the universe is "created from noth- 
ing" (Zel'dovich k Grishchuk 1984). The universe 
which emerges from nothing, created by a quantum 
tunneling process, is generically quite homogenous if 
the tunneling probability is suppressed enough (which 
selects the most symmetric path as the most proba- 
ble one). The tunneling thus solves the homogeneity, 
isotropy and horizon problems before "inflation" be- 
gins, allowing for a truncated period of inflation so 
as not to inflate the curvature radius well outside the 
horizon. Specific models which lead to closed uni- 
verses can be constructed (Linde 1995) using scalar 
field (inflaton) potentials of the form <f) 2 exp(c/> 2 ) or 
by introducing non-minimal couplings of the inflaton 
field to the curvature R. We shall assume for definite- 
ness that inhomogeneities in the bubble which forms 
our present universe by tunneling are smaller than 
those produced during the subsequent period of in- 
flation (Linde k Mezhlumian 1995), though this need 
not necessarily be the case. 
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Figure 1: Contours of the age of the universe, in units 
of the Hubble time Hq 1 . Also shown (solid line top 
left) is the critical line for an infinite age universe and 
(solid line bottom right), a universe which becomes 
asymptotically static. Universes above (below) these 
lines have no big bang (have a big crunch). Lines 
of constant qo = \£l m + £l r — are parallel to the 
q = line shown dotted. 



4 The Classical Tests 

A thorough discussion of the classical cosmological 
tests in a closed universe can be found in e.g. Sandage 
(1961,1988), Gott, Park k Lee (1989), Durrer k 
Kovner (1990), Bjornsson k Gudmundsson (1995). In 
this section we remind the reader of some properties 
of closed universes which will be necessary to under- 
stand our later results. Helpful diagrams of a(t) are 
presented by Felten k Isaacman (1986). Readers who 
are unfashionable enough to be familiar with closed 
models can skip this section entirely. 

Except in §8, where we consider an extreme model, 
we will restrict ourselves to the parameter range 
(Q m ,Q\) G [0,2] x [0,2]. We will use the notation 
fi A = A/3.H", 2 for the effect ive contribution of the 
cosmological constant to fi. The general Friedmann 
equation is 

— = Ho\jQ, r a~ i + Q, m a~ 3 + SJifd -2 + Q,\ 
a 

= H E(z) , (6) 
with 1 + z = a -1 and oq = 1. The acceleration equa- 



tion can be written 

^ = -H 2 ^a- 3 + n r a- 4 -n A ). (7) 

Since fi r is so small today, we make little error by 
neglecting it for this discussion, although we always 
include it in numerical calculations. 
Using (6), the age integral is 

f°° dz' 

H t(z)= / — — — , (8) 

1 } J z (l + z>)E(z>y 1 } 

while the conformal time integral is 




We show contours of ifo^o in Fig. 1. As can be 
seen, increasing £l\ for fixed fi ra results in longer 
ages. This is one of the reasons why cosmological con- 
stant models are re-considered periodically (e.g. Gunn 
k Tinsley 1975, Tinsley 1977, Zel'dovich k Sun- 
yaev, Sandage k Tammann 1984). While the age (in 
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Hubble units) increases rapidly towards the top left of 
Fig. 1, we will see that much of that region of param- 
eter space is ruled out by classical cosmological tests 
and by lensing constraints. For a concrete feeling of 
what these ages mean, we can write 



12.22 Gyr (H t ) 



Ho 
80 



(10) 



By considering the acceleration equation (a = 0), 
we see that there is an inflexion point in the Hubble 
constant evolution when 



(1 + Zloitf 



(11) 



This 'loitering' phase can be tuned to be arbitrarily 
long, giving rise for example to the large ages men- 
tioned above. For models with even modest amounts 
of loitering one runs into problems with classical cos- 
mological tests which measure volumes. In the limit 
of extreme loitering, one even "recovers" Olber's para- 
dox. 

The age of the universe becomes infinite in a model 
where a — > at the inflectional point, or equivalently 
where E(z) is zero at its extremum. These two con- 
ditions together imply 



27fi^fi A = 
This has solution 

riyv = i — 



4(n n 



I) 3 - 



(12) 



{ 



-fi 2 / 3 x 
2 m 



1/3 



1 — fi ra — \J\ — 2fi r 

1 - n m + y/\ - 2fi ra ) J , (13) 



which can also be written in other forms (see Carroll, 
Press & Turner 1992). We show this "critical line" in 
the top left of Fig. 1. 

These infinite age universes are also known as 
Eddington-Lemaitre models. The degenerate case 
f2 ra = 0, £l\ = 1 is the de Sitter or Steady State 
model. Models near to this case have a period in 
the past where a ~ 0, and are therefore known as 
the "loitering", "coasting" or "hesitating" models. In 
the past they have also been called Lemaitre mod- 
els, and have long been popular as a means of im- 
proving the ability for galaxies to have formed by the 
present (Lemaitre 1931, Rawson-Harris 1969, Brecher 
k Silk 1969, Tomita 1969, Heath 1977, Occhionero et 
al. 1980); we shall discuss the issue of growth in loiter- 
ing models further in §6. A large amount of effort was 
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Figure 2: Solid lines show contours of lensing prob- 
ability, relative to the Einstein-de Sitter model for a 
source at redshift 2. Dashed lines show the antipode 
redshift. Both P\ ens (z s = 2) > 5 (probably a con- 
servative limit, although estimates vary over ~ 2-10) 
and z a < 4 are ruled out. The lensing constraint is 
clearly the stronger. 



devoted to these Lemaitre models in the early days of 
high redshift quasars and cosmological number counts 
(e.g. Shklovskii 1966, Kardashev 1966, McVittie & 
Stabell 1967, Petrosian, Salpeter & Szekeres 1967, 
Solheim 1968, Petrosian 1969, Petrosian 1974). To- 
day the lack of a pile-up at any redshift rules out 
our being too close to the Eddington-Lemaitre case. 
Note that Loitering universes have divergent r\ as well 
as t. In fact any "path length" integral diverges, since 
E(z) — > for an arbitrarily long range in z (e.g. the 
lensing probability shown in Fig. 2). 

Models with higher £l\ than the critical line have 
no Big Bang, i.e. they contract from infinity, bounce, 
and then go into a period of exponential expansion. 
Such models therefore have a minimum value of a, 
which means a maximum observable redshift. In fact 
such models must have 



3) 



(14) 



(Crilly 1968), and since we observe objects out to z ~ 
5, while fi ra 0.01, such models are not allowed 
(Ehlers & Rindler 1989). Hence the upper left corner 
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of Fig. 1 is definitively ruled out. 

Models which have £l\ = and fi ra > 1 will ulti- 
mately collapse. However, because of the scale-factor 
dependence in the Friedmann equation, most models 
with £l\ > will become dominated by the cosmolog- 
ical constant and end up in a de Sitter phase, i.e. they 
will have an eternally inflating future. There is a small 
region of parameter space for which fi ra wins, and the 
universe collapses before A can take over. The criti- 
cal line is when a — > for the (future) inflection point 
a = 0. Hence this line is also a solution to (12). 
Specifically the critical line for collapse is 




for fi ra > 1 (Felten k Isaacman 1986). This is shown 
by the solid line in the lower right of Fig. 1. 
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5 The CMB 

The postitive curvature of a closed model can shift 
features in the CMB fluctuation power spectrum to 
larger angular scales. It has been conjectured for some 
time (Blanchard 1984, Durrer k Kovner 1990) that 
this might lead to a reduction in the amplitude of the 
anisotropics, however no explicit calculations of the 
anisotropy spectrum have been done. 

Early work on CMB anisotropics in non-flat geome- 
tries (e.g. Wilson 1983, Tomita k Tanabe 1983, Ab- 
bott k Schaefer 1986) concentrated on the large-angle 
part of the spectrum. In this section we present the 
first full calculations of the angular power spectrum 
of fluctuations in a closed universe. The formalism 
employed is reviewed in the appendix, where the Ein- 
stein, fluid and Boltzmann equations for an arbitrary 
Friedmann-Robertson- Walker geometry are derived. 
The results were obtained by numerical solution of 
these coupled equations (see HSSW, and references 
therein, for more details). 

Conventionally one expands the temperatures on 
the sky in terms of a complete set of functions: 

T{p,<t>) = T CM B^a lm Y lm {p,<t>), (16) 

lm 

where the spherical harmonics Yi m are the curved sky 
equivalent of Fourier modes, and Tcmb is the average 
temperature of the CMB, so that the a£ ra 's are di- 
mensionless. The power spectrum is then defined in 
terms of the ensemble average of the coefficients: 

Q = (K ra | 2 ). (17) 



Figure 3: Contours showing the location, i'peak ) of the 
first peak in the CMB anisotopy power spectrum. The 
other peaks in the power spectrum form an almost 
harmonic series at 2^ pe ak> 3£p e &k--- Note the region 
in the upper right of the plot is allowed by classical 
tests, but will be ruled out by the CMB if the peak is 
at i'peak ~ 200 as suggested by recent data. 

The power per logarithmic interval in £ is approxi- 
mately 1(1 + l)Ci, and this vs \og£ is what we shall 
refer to henceforth as the angular power spectrum. 
For reference, multipole £ probes an angular scale of 
~ 60/£ degrees (see appendix of White et al. 1994) 
for all but the smallest i"s. 

The most obvious feature of the angular power 
spectrum is the peak that occurs around £ ~ 200 (see 
Fig. 4). This peak represents an acoustic oscillation 
that was maximally overdense at last scattering (for 
a discussion see Scott k White 1995, Hu, Sugiyama 
k Silk 1995). Recall that for fi tot ~ 1 and z > 1 the 
angle subtended by the horizon is roughly (fitot/ z) 1 ! 2 
(e.g. Weinberg 1972, eq. [15.5.39]). Since features in 
the angular spectrum of CMB fluctuations are related 
to the horizon-size at z ~ 1100, we expect that the 

— 1/2 

whole spectrum will shift oc f2 tot . So for example, 
the main acoustic peak of adiabatic models will be at 

^peak — 220 ^tot 

One can estimate this effect more precisely by cal- 
culating the angular scale subtended by the "sound 
horizon" at last scattering. The first peak appears at 
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(Hu & Sugiyama 1995) 



<-peak 



(18) 



where, using xo — X* = VK (r]o — »7*), we define 

r* = sin [a/k^o - Ilss)] , (19) 

the effective distance to the last-scattering surface. 
The sound horizon is 



r s (v) 



dr] c s 



(20) 



o 



where c s is the speed of sound. An analytic expres- 
sion for r s can be found in Hu k Sugiyama (1995). 
The higher peaks occur in an almost harmonic series, 
i.e. 2i ? p ea j i , Si^pgak, . . . 

For the case of an open universe, the sine function 
of (19) is replaced by sinh, which is monotonic in its 
argument. For closed models, however, there is the 
possibility of the argument getting bigger than 7r/2, 
or of being mr (i.e. an antipode) in which case r* can 
have more interesting behavior. 

We show in Fig. 3 the position of the first peak as a 
function of fi ra and Here we have assumed that 
recombination occured at z rec = 1100 for all mod- 
els, since there is little variation with cosmological 
parameters. In the numerical work we have followed 
the recombination process accurately (HSSW). 

Notice that towards the top left corner of Fig. 3 the 
peak moves from high to low £ and back repeatedly. 
This is a consequence of the last scattering redshift se- 
lecting the equator, the antipode, the equator and so 
on. In the familiar A = closed models (5), a light ray 
leaving the Big Bang arrives back at the same point 
at exactly the Big Crunch (x = 2ir). However, this 
is no longer true for A > closed models, where it is 
possible to have x ^ 27r, i.e. many antipodes rather 
than just one. [This is easy to see for an extreme 
loitering model, where the universe remains static 
for a long period, allowing light rays to circumnav- 
igate the universe many times.] The peculiar cosmo- 
logical consequences of an observable antipode have 
been investigated by many authors (e.g. Petrosian k 
Salpeter 1968, Rowan- Robinson 1968, Baylis 1970, Bi- 
raud k Mavrides 1980, Gott 1985, Gott k Rees 1987, 
Gott, Park k Lee 1989, Durrer k Kovner 1990, 
Bjornsson k Gudmundsson 1995), although the pos- 
sibility now seems quite remote. The redshift of the 
first antipode, z a , was shown in Fig. 2. 

The angular power spectra for a variety of models 
are shown in Figs. 4-6. One can see that the posi- 
tion of the peak moves according to the prescription 
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Figure 4: Angular power spectra of CMB anisotro- 
pics for models with fi ra > 1. The position of the 
first peak moves right with decreasing fi ra . We have 
chosen h = 1 and fi^/i 2 = 0.0125 for simplicity. Note 
the steep rise of the spectrum towards low £, which 
allows us to rule out high models. 

outlined above. For these models we have specifically 
assumed that h = 1 and = 0.0125. The relative 
heights of the peaks will depend on this choice (higher 
£Ib will enhance the first peak for example) but the 
positions of the peaks and the large angle (small €) 
anisotropy spectrum will not. 

In Figs. 4-6, we have assumed that the primordial 
fluctuation spectrum is "scale invariant", or that the 
potential fluctuations, <J>, are constant per logarith- 
mic interval in wavenumber (see appendix). Such a 
spectrum is more physically motivated (e.g. from in- 
flation) than assuming the fluctuations are a power- 
law in wavenumber, since in a non-flat geometry there 
is an important scale in the problem: the curvature 
scale. Moreover, such a primordial spectrum leads 
to smaller deviations from flat CMB spectra than 
the corresponding power-law. We show the primor- 
dial spectra vs. wavenumber in Fig. 7, where one can 
see that the large scale power is reduced (over the 
power-law case) in a closed model and enhanced in 
an open model, counteracting the geometric tenden- 
cies to have more (less) large angle power in closed 
(open) models. The apparent divergence of P pr i m at 
small wavenumbers in the open model is counteracted 
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Figure 5: Angular power spectra for models with 
fltot = 2 and £l\ > 0. These models rise more 
steeply towards low £ than the £l\ = curves, show- 
ing that for constraining fltot> the most conservative 
constraints are with £l\ = 0. 

by the volume element, leading of course to finite fluc- 
tuations. 

In the top right of Fig. 3 the peak location is 
^peak = 0, which is just where the last scattering sur- 
face lies at the antipode. It had been conjectured 
that this would be a way to suppress CMB anisotro- 
pics (Durrer & Kovner 1990). We show in Fig. 6, 
the angular power spectrum taking into account the 
finite duration of last scattering, and the non-trivial 
gravitational effects on the spectrum, for a model with 
^peak = according to (18). Notice that the spectrum 
is featureless, showing only the "damping" tail of the 
anisotropics. However, it is non-zero because a fi- 
nite width last scattering surface cannot project onto 
quite the whole sky, and furthermore there are still 
gravitational interactions on the photons. So we still 
see some anisotropics, but the acoustic peaks have 
all been magnified to ultra-large angular scales (very 
small £). 

We also show on Fig. 6 two models where the last 
scattering epoch occurs slightly before and slightly 
after the antipodal redshift. Note that as we shift 
the antipode through z = 1100, the structure in the 
spectrum moves to £ ~ and back out to finite £ 
again. For all of these models the predicted large an- 
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Figure 6: The solid line shows the angular power spec- 
trum for a model with antipodal redshift at the last 
scattering surface: z a = z rec (specifically: fi ra = 0.5; 

= 1-43). Note the absence of any structure in the 
spectrum. With z a = z rec all structure is projected to 
£ = 0, and we see only the damping tail plus some in- 
tegrated Sachs- Wolfe effects at higher £. Also shown 
are two models where the peak is moved towards £ = 
(Aa = 1-35; dotted) and where it has moved through 
£ = back to higher £ (Q\ = 1.5; dashed). 

gle anisotropy (e.g. C2) for fixed large-scale matter 
fluctuations is within a factor of 10 (in power, so ^/To 
in temperature) of the standard cold dark matter re- 
sult. Thus, though the power falls off rapidly with 
£, the suppression in power is not as extreme as one 
could naively imagine. 

From Figs. 4-5 we see that on large scales (small 
£) the spectra are steeply falling functions of £. This 
was first noted by Abbott & Schaefer (1986), and is 
due to the decay of the gravitational potential at late 
times. A photon falling into a potential well which 
is decaying with time, has its in-fall blueshift only 
partly cancelled by its climb-out redshift. The net 
anisotropy builds up along the path, leading to larger 
temperature fluctuations on large scales (or long char- 
acteristic times where the effect has the most time 
to operate). This is known as the (late) integrated 
Sachs- Wolfe effect (Sachs & Wolfe 1967; see White et 
al. 1994, Hu, Sugiyama & Silk 1995, Hu & White 1995 
for more discussion). Because this is the dominant 
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Figure 7: The primordial power spectrum for con- 
stant gravitational potential fluctuations, <J>, per log- 
arithmic interval in wavenumber, also called a "scale 
invariant" spectrum. The points on the fi ra = 1.5 
curve indicate the quantized values of wavenumber, 
/3 = 3, 4, 5 . . ., in a closed universe. Note that only 
the fi ra = 1 model looks like a power law when plot- 
ted in this way, and that the behaviour for fi ra > 1 
continues from the fi ra < 1 region as expected. 

effect at large scales, and does not appear in a uni- 
verse with critical matter density, it makes scaling ar- 
guments difficult without accurate calculation of the 
anisotropics, evolved to the present epoch (c.f. Bunn, 
Scott k White 1995, White k Bunn 1995, Sugiyama 
1995, Stompor, Gorski k Banday 1995). 

As can be seen from Fig. 3, the location of the 
first peak in the power spectrum can allow one to 
probe a large region of the ^l m -^l\ plane. Data on 
anisotropics on degree scales are not yet up to the 
task of strongly constraining this position (Scott, Silk 
k White 1995), however we can still place new con- 
straints on fltot from the large-angle fluctuations mea- 
sured by the COBE satellite. We shall return to this 
in §7. 

If however you believe in the existence of an acous- 
tic peak in the anisotropy data, then this does place a 
strong constraint on closed models. It seems unlikely 
that the current data can be fit without some sort of 
peak at £ ~ 100-200, unless several experiments have 
not been measuring primordial anisotropics. Taking 



the optimistic view that the data will indicate this 
more definitively in the near future, we may soon have 
a constraint that, for example, fpeak > 100. From 
Fig. 3 we can see that this will rule out a substantial 
region of parameter space. We shall see in §7 that 
this region is probably already ruled out by COBE 
constraints alone. However, a tight constraint from 
the peak position would be insensitive to details of 
the primordial power spectrum around the curvature 
scale, where it is perhaps most dependent on specific 
inflationary models. Eventually we would expect an 
accurate measurement of the anisotropy spectrum to 
yield an extremely precise determination of fltot- 

6 Relation to Large Scale Structure 

In Fig. 8 we show the relative normalization of the 
large-scale matter power spectrum and large-angle 
CMB angular power spectrum for a variety of models. 
Since much of the large angle CMB spectrum comes 
from the integrated Sachs- Wolfe effect, this normal- 
ization ratio cannot be obtained from simple growth 
factor arguments alone. We have calculated numer- 
ically the ratio of B to C±o where C±o is as above, 
while 

P(P) 



B 



prim 



(/3) 



(21) 



/3=3 



is the primordial power spectrum amplitude. We have 
defined this as the analogue of the quantity B some- 
times used to normalize the CDM power spectrum for 
fi ra = 1 (e.g. Bunn, Scott k White 1995). There B 
is defined through P{k) = BkT 2 (k) (for an n = 1 
spectrum) where T(k) is the transfer function. Using 
T(k — > 0) = 1 we can write B = limj^o P(k)/k, 
from which (21) follows naturally, once you notice 
that a closed model has wavenumber fiyfK quantized 
in units of the curvature scale (see appendix for de- 
tails). Since the /3 = 1 and /3 = 2 modes are pure 
gauge modes (Abbott k Schaefer 1986) the smallest 
wavenumber (largest physical scale) is /3 = 3. We 
have chosen to use Cio, rather than a lower £, to min- 
imize the primordial power spectrum dependence and 
reduce the contribution from the integrated Sachs- 
Wolfe effect. For £ = 10 we are still very insensitive 
to details of cosmological processing and parameters 
other than fi. 

For the matter power spectrum we also focus on 
the largest scales, to avoid complications due to the 
variable redshift of matter-radiation equality, which 
affects the amplitude of the fluctuations on smaller 
scales. As Fig. 8 shows, the matter-radiation normal- 
ization is a smooth function of fl t ot which depends 
upon the growth of fluctuations between z = 1100 
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Figure 8: The relative normalization of the large- 
angle CMB and large-scale matter power spectra in 
closed models (see text for details). B is a mea- 
sure of the large-scale normalization of the matter 
power spectrum, while Cio measures CMB fluctua- 
tions on scales smaller than the curvature scale but 
large enough to probe mostly potential fluctuations. 

and the present. For those cases we considered, which 
cover some extreme universes, the ratio of the mat- 
ter to radiation power spectra varies by a factor of 
~ 10 (cf. Sahni, Feldman & Stebbins 1992, and many 
earlier discussions of exponential growth during a loi- 
tering phase). The ratio for £l\ = and 1 < fi ra < 2 
is well fit by 

= 320 exp (-1.16 n m ), (22) 

Cio 

in units of (10 5 h~ x Mpc) 3 . The relation between the 
CMB and the smaZZ-scale matter fluctuations is given 
by Fig. 8 and the transfer function T(k), to which we 
now turn. 

For flat and open models it is known that the trans- 
fer function, T(k), relating the initial and final power 
spectra is a function simply of the redshift of matter- 
radiation equality. Often one uses the fitting function 
(Efstathiou 1990) 

T(k) = [l + {ak + (bk) 3 ' 2 + (ck) 2 Yj 1 , (23) 

with parameters a = (6.4/r)Mpc, b = (3.0/r)Mpc, 
c = (1.7/r) Mpc, v = 1.13 and T ~ tt m h. Notice that 



Figure 9: The matter transfer function. The solid 
lines are from the fitting function described in the 
text, while the symbols show the results of explicit 
calculations. Apart from a small discrepancy near 
the turn-over, which is caused mostly by the baryon 
content assumed for our models, the curves are a good 
fit to the numerical calculation. 

perturbations on small scales, specifically those which 
entered the horizon when the universe was radiation 
dominated, have had their growth suppressed relative 
to those on larger scales. 

We show in Fig. 9 that (23) is a good fit to the 
closed models as well. Thus, even though the growth 
rate can be "complicated" by the presence of loitering 
phases, T(k) today is still fit by a 1 parameter fam- 
ily of curves. It is true that there are modes which 
can grow quickly during the loitering phase, when 
matter and cosmological constant conspire to keep 
the universe almost static as in the original analysis 
of (Jeans 1928). However, they will also have their 
growth retarded when the cosmological constant be- 
comes dominant at later times and the universe be- 
gins to expand exponentially. These two effects ap- 
proximately cancel, leading to fluctuation amplitudes 
vs. wavenumber that would be predicted from analogy 
with the open or flat models. [A similar conclusion 
for the large-scale (small k) modes was reached by 
Durrer k Kovner (1990).] 

Using the large-scale normalizations of Fig. 8 and 
(23) one can calculate the moments of the power spec- 
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tram for whatever parameters (fl ra , h) are of in- 
terest. The "window functions" for moments such 
as as, V Tms , differ from those in a flat universe (see 
e.g. Wilson 1983) by geometric factors related to the 
3- volume of the metric (i.e. r — > sinx)- 

7 Limits on from the CMB 

There are already some weak constraints on flat A- 
dominated models from the shape of the large-angle 
power spectrum, amounting to £l\ < 0.8-0.9 (Bunn 
& Sugiyama 1995, Stompor, Gorski & Banday 1995, 
White & Bunn 1995). We can similarly use the 2-year 
COBE data to put a limit on fl t ot from the shape of 
the CV's at small £. This is quite insensitive to the 
cosmological parameters other than fltot> making the 
limit reasonably robust. 

First note from Fig. 5 that for fixed fl t ot the spec- 
tra with > fall more steeply than with £l\ = 0. 
This is because there will always be more ISW ef- 
fect in a A model (Kofman & Starobinsky 1985), at 
redshifts where the cosmological constant starts to 
dominate: z\ ~ ^A^tot — I) -1 - Since the COBE 
data prefer spectra which are approximately "flat", 
or C~[ x oc 1(1 + 1), it follows that for fixed fl t ot mod- 
els with Oa > will be less likely than the model 
with Oa = 0. Also models with initial power spec- 
tra which have constant power in gravitational po- 
tential per logarithmic interval in wavenumber (scale 
invariant spectra) are "flatter" than models whose 
power spectra are proportional to wavenumber (see 
appendix and Fig. 7), and thus will be preferred by 
the COBE data. 

In Fig. 10 we show the likelihood function vs. fl ra , 
integrated or marginalized over the normalization, for 
models with scale invariant spectra and Oa = 0, fit 
to the COBE data. The likelihood function is nor- 
malized so that a "flat" spectrum, C'^ 1 oc 1(1 + 1), 
has £ = 1, and we have included a "prior" that 
fi ra > 0.1 (i.e. we set £ = for fi ra < 0.1). The 
likelihood function for fl ra < 1 comes from (Gorski et 
al. 1995) rescaled to the normalization of (White & 
Bunn 1995). The dashed line shows the cumulative 
likelihood, from which we estimate the 95% CL up- 
per limit on fl ra to be fl ra < 1.5. Any model which 
has fltot > 1-5 will be less likely than the model with 
fl ra = 1.5 and Oa = 0. If we take the upper limit 
from the cumulative likelihood to mean fi ra = 1.5 is 
statistically inadmissable at the 95% CL, we obtain 
the limit fl t ot < 1-5. 

This limit on fl t ot compares well with other limits 
on this quantity (Loh & Spillar 1986), but it is de- 
pendent on the initial power spectrum assumed. If 
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Figure 10: The likelihood function (solid) vs. fl ra , for 
models with scale invariant primordial spectra and 
Oa = 0, fit to the COBE 2-year data and marginalized 
with respect to normalization. We have set £(£l m ) = 
for fl ra < 0.1 (our "prior"). The dashed line shows 
the cumulative likelihood, and the arrow indicates the 
95% CL upper limit on fl ra . As discussed in the text 
this is a conservative upper limit on fltot- The results 
for fl ra < 1 come from Gorski et al. (1995). 

the primordial power spectrum is "redder" than scale 
invariant (e.g. n < 1) then our limit is strengthened. 
However should a particular inflationary model pre- 
dict n > 1 our limit would be weakened. A less model 
dependent limit could be placed on fl t ot from the po- 
sition of the peaks in the angular power spectrum, 
when the data becomes less ambiguous. 

8 An Extreme Model 

As an example of how non-standard models can 
sometimes be a surprisingly good fit to data, (Har- 
rison 1993) considered a model with fi ra = 10 and 
H = 10kms -1 Mpc -1 . The idea was that the Hub- 
ble flow in such a universe could be so chaotic that 
our local estimates of Ho could be far from the global 
value. With h = 0.1 the model is very old, while 
with fi ra = 10 it is very small and has almost finished 
its expansion and begun its (re)collapse. While we 
consider this model to be far from acceptable when 
examined in detail, we will use it here as an example 
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Figure 11: The angular power spectrum for the model 
of Harrison (1993) with fi ra = 10 and h = 0.1 as- 
suming a scale invariant primordial spectrum and 
tt B h 2 = 0.0125, tt A = 0. Notice the peak in the 
power spectrum is at £ ~ 100, while the large-angle 
structure of the power spectrum is wildly inconsistent 
with the COBE data. 

of an extreme closed model. Harrison presents some 
scaling arguments for how structure formation and 
microwave background anisotropics would work out 
in such a model. Here we calculate P(k) and the Ct 
spectrum in detail. 

We shall assume that = 1.25 so that the model 
agrees roughly with Big Bang Nucleosynthesis deter- 
minations of fifi/i 2 , though we can relax this assump- 
tion without changing our principle results. Since we 
know observationally that AT/T ~ 10 -5 on large 
scales (Smoot et al. 1992) we can have confidence that 
our calculation of the angular power spectrum in lin- 
ear perturbation theory will be accurate even in this 
extreme model. 

The resulting anisotropy spectrum is presented in 
Fig. 11 for a scale invariant primordial power spec- 
trum. One can see that the first peak occurs at 
£ ~ 100 (cf. Fig. 3; the height of the peak will depend 
on our assumption of the baryon content, but not its 
position) and at large angles it is wildly inconsistent 
with COBE. In linear theory the large-scale density 
fluctuation to large angle CMB fluctuation ratio is 
B/C 10 ~ 10 (10 5 h' 1 Mpc) 3 . Thus, normalized to 



COBE, the large scale power in this model will be 
less than in e.g. standard CDM. 

Once again the linear theory transfer function from 
explicit evolution of the density perturbations agrees 
well with (23). We find slightly less small scale 
power in the numerical calculation than for (23) with 
T ~ £l m h = 1, but only at the ~ 10% level for 
k < 0.2/i/Mpc. A model with T ~ 1 has much more 
small scale power than a standard CDM model, lead- 
ing to earlier non-linearity at small scales. It would 
probably be in conflict with existing measurements of 
the large-scale matter power spectrum, which appear 
to require T ~ 0.25 (see e.g. Peacock & Dodds 1994). 
If we had decided to have = fi ra = 10 (in vi- 
olation of nucleosynthesis bounds on £Ib) then the 
transfer function would damp much more strongly at 
high-fc than for the model with 90% dark matter, and 
would show oscillations similar to those which appear 
in the CMB anisotropy power spectrum. 

9 Other Limits 

The observational limits on £l\ have been reviewed 
by (Carroll, Press & Turner 1992), so in this section 
we simply make a few additional remarks, referring 
the reader to that work for more details. Existing 
and future limits that are comparable to those from 
the CMB are discussed below. 

Galaxy counts once seemed like a good way of mea- 
suring £l\ through the change in the volume element 
(e.g. Loh & Spillar 1986), but separating the effects of 
geometry and evolution has proven to be the downfall 
of such methods (Gardner, Cowie & Wainscoat 1993). 
Nevertheless, the results from number counts have 
been interpreted as implying that fi ra < 1.5 (Loh & 
Spillar 1986), which is similar to our limit. 

As pointed out in Gott, Park & Lee (1989), the 
existence of "normal looking" quasars at z > 4 
puts a lower limit on the redshift of the antipode: 
^/K[t]o — *?( z a)] = I"- The contours of antipodal red- 
shift are shown in Fig. 2, where we see that z a > 4 
rules out a small region in the top left of the plot. As 
discused in §5 the presence of structure in the CMB 
anisotropy spectrum would put a much stronger limit 
on the redshift of the antipode: z a > 1100. The pos- 
sible loop hole in this argument, namely that the last 
scattering surface lies beyond the antipode so that 
the structure "reappears" at £ > 0, will be seen to be 
ruled out by gravitational lensing constraints. How- 
ever we will see that these constraints are in any case 
stronger than the lower limits on z a . 

In a universe which becomes £l\ dominated at some 
redshift, z\, the expansion rate is dramatically differ- 



12 



ent before and after z\. Such a change can manifest 
itself as a "break" in the number of Ly a absorbers 
per unit redshift (Fukugita k Lahav 1991, Turner k 
Ikeuchi 1992). Future surveys of low- z Ly a clouds 
may be able to put strong constraints on (or prefer) 
models, but at present the situation is unclear. 
Another sensitive probe of £l\ in the future may 
come from separations of close quasar pairs (Phillipps 
1994), as measured for example through the clustering 
properties of quasars in redshift surveys such as the 
Sloane Digital Sky Survey or the AAT 2dF. Here the 
idea is that the underlying correlation properties of 
quasars are expected to be isotropic, but when mea- 
sured in terms of angular position and redshift one 
expects distortions (e.g. in the familiar tt — a plane 
for various redshift shells). The comoving separation 
of a pair of objects at z±6z and separated on the sky 
by 60 is (Phillipps 1994) 



g\8zf 



J2 (zq^ 

a » 



l M 
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(24) 



where cIm(z) is the "proper motion" distance (Wein- 
berg 1972) and 



g(z) = (1 - [n tot - l]H 2 d 2 



-1/2 (IcIm 



Ml 



(25) 



The ratio of cIm to its derivative, g, at high redshift, 
which is a measure of the distortion of the isotropic 
correlation function, is strongly dependant on the cos- 
mology assumed. With 10 5 quasars this method could 
result in extremely strong constraints on (or a detec- 
tion of) A, but awaits data from large surveys. 

The strongest limits on £l\ to date come from the 
rarity of gravitational lensing events (e.g. Fukugita, 
Futamase k Kasai 1990, Turner 1990, Fukugita et 
al. 1992, Kochanek 1993, Carroll, Press k Turner 1992). 
In a universe with a loitering phase the path length to 
any given redshift is dramatically enhanced, and the 
expected number of gravitational lenses along any line 
of sight (assuming a fixed density of objects) far ex- 
ceeds the observed number. Define P\ ens (z s ) as the 
probability of lensing a source at redshift z s relative 
to that in an fi ra = 1 universe (Carroll, Press k 
Turner 1992, Fukugita et al. 1992). The current sur- 
veys require P\ ens (z s ~ 2) < 2-10 depending on the 
detailed assumptions and data used; we shall take 
-Piens < 5 to be a conservative limit. Contours of 
-Piens(^s = 2) are shown in Fig. 2. 

Since the lensing limit is currently more stringent 
than all of the other limits, the allowed region of pa- 
rameter space becomes that to the lower right of some 



2-10 line, and lower left of fit 
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model dependent bounds on fl t ot could come from the 



Figure 12: Constaints in the ^l m -^l\ plane. The un- 
shaded region is currently allowed by fairly conserva- 
tive constraints. The three sides of the region come 
from lower limits on fl ra , upper limits on Pi en s and 
an upper limit on fl t ot from large-angle CMB aniso- 
tropics. We also show two age lines with -ffo^o = 2/3 
and Hoto = 1. Note that a fairly wide range of closed 
models have not yet been definitively ruled out, and 
that some of them are relatively old. 



location of the peak in the anisotropy spectrum, once 
the degree scale data becomes more robust. Taking 
-Piens < 5 it is still possible to have Hoto ~ 1 for a 
small region of parameter space with fl ra < 0.5. The 
lensing and CMB constraints are acting in different 
directions, thus an "improvement" in either will help 
to squeeze the allowed region in parameter space. 

We show the currently allowed region of parameter 
space in Fig. 12, using dynamical, lensing and CMB 
constraints. We have chosen to be reasonably conser- 
vative in each case. The allowed region still contains 
a variety of closed models. We have indicated two 
representative age lines on this plot also. Note that 
for Ho = 80 kms _1 Mpc _1 , the Hoto = 1 line has an 
age of 12.2 Gyr, while the Hoto = 2/3 line has an age 
of only 8.2 Gyr. Since dynamical estimates tend to 
give values like fio ~ 0.3, it is perhaps tempting to 
speculate that the best bet for a model which is old 
enough, may be one with fl ra + £l\ > 1. However, it 
is also clear from the figure that the age gain over a 
fltot = 1 model is not very great. 
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10 Conclusions 




Figure 13: The position of the first peak in the CMB 
power spectrum, ^ pe ak> as a function of fi ra for -ffo^o = 
1 (solid), 1.1 (dotted), 0.9 (dashed), fpeak is a rapidly 
varying function of fi ra along lines of constant .Hero- 
in the future we expect strong constraints on -ffo^o 
from a combination of large-scale structure (measur- 
ing fl ra ) and the CMB. We show an example of this 
in Fig. 13, where we have plotted the position of the 
first peak in the CMB power spectrum as a function 
of fl ra , for three values of -ffo^o- Note that the posi- 
tion of the peak in the experimentally favoured region 
is a rapidly varying function of fi ra . 

In closing we note that large scale dynamics are 
mostly insensitive to A (see e.g. Martel & Wasser- 
man 1990, Martel 1991, Lahav et al. 1991). So, for 
example, the comparison of density and velocity fields 
measures the function fi^ 6 /6 (where b is the "bias" 
or amplitude of the matter power spectrum) indepen- 
dently of whether the model contains a contribution 
from A which makes it flat. In the limit that the de- 
pendence of the perturbation growth rate is a weak 
function of these conclusions apply equally well 
in closed models. Dynamical tests thus generally tend 
to measure fi ra rather than, say, fltot- I n marked con- 
trast the location of the peaks in the CMB anisotropy 
angular power spectrum depends essentially only on 
the geometry, and thus measures both fi ra and 
largely through the quantity fltot- This makes the 
combination of the two an extremely powerful con- 
straint on cosmological models. 



We have examined CMB anisotropics and large- 
scale structure formation in models with positive spa- 
tial curvature, i.e. closed models. We presented for 
the first time the angular power spectrum of the an- 
isotropics over a range of angular scales. The relative 
normalization of the CMB anisotropics and the large- 
scale matter power spectrum show behaviour which 
is a simple continuation from the flat case. 

For scale-invariant potential fluctuations, the an- 
isotropics are steeply rising towards low multipoles. 
With A = the position of the main acoustic peak 
scales approximately as fpeak cx fi to {/ 2 as expected. 
There is a richer structure (shown in Fig. 3) when 
A ^ 0. It is possible in principle to have last- 
scattering occur at the antipode, in which case the 
spectrum looks very different. However, this case is 
unlikely to occur for our observed universe. 

As a particular example we have calculated the 
anisotropy spectrum for the model postulated by Har- 
rison (1993), which has fl ra = 10 and h = 0.1. Since 
the CMB fluctuations are measured to be well in the 
linear regime, the chaotic behaviour envisioned by 
Harrison to account for our "misinterpretation" of the 
current data is hardly relevant, so one can say with 
confidence that the model is in conflict with current 
CMB and large-scale structure data. 

Extending the likelihood function vs. from the 
COBE 2-year data to > 1 we obtain a 95% CL 
upper limit of fl t ot < 1-5 for primordial spectra with 
n < 1, which compares favorably with existing limits 
on fi. This limit is dependent on our assumption 
of a primordial power spectrum and the absence of 
tensor modes (which are self-consistent), however a 
more secure bound could come from the position of 
the first peak in the angular power spectrum as shown 
in Fig. 3. 

The low number of gravitationally lensed quasars is 
currently the tightest constraint on the value of 
However, limits from large-angle CMB anisotropics, 
as well as the position of the acoustic peak, show 
great promise for the near future. It is also impor- 
tant to note that these effects can constrain different 
combinations of fi ra and £l\ than lensing and large- 
scale structure (which measures largely fl ra ). If we 
were able to measure fi ra accurately (from the mat- 
ter power spectrum for example), then for fi ra < 0.5 
an "old" universe has a strong dependence of fpeak on 
fi ra (see Fig. 13). Thus it should be easy to measure 
either £l\ or -ffo^o from the peak location in such mod- 
els, given the other. This is just one consequence of 
the age contour lines running roughly perpendicular 
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to i'peak contour lines in the (fl ra , plane. 

There is still a small region of parameter space 
which is allowed by cosmological constraints, and 
which lies in the fl t ot > 1 part which is oftenneglected 
(Fig. 12). For some smaller part of this region the age 
of the universe is also favourable for the high values of 
the Hubble constant which seem most commonly to 
have been discussed in the last few years. With this 
motivation, but largely because they are still viable 
today, we believe these closed models are worthy of 
examination. 

We wish to thank Wayne Hu for many useful con- 
versations and Ted Bunn for help with the COBE 
data. 
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A Perturbations in Non-Flat Geometries 



In this appendix we provide a derivation of the Einstein, fluid and Boltzmann equations in a Friedmann- 
Robertson- Walker universe with arbitrary spatial curvature, paying particular attention to the case of a closed 
universe (spherical geometry). Our derivation is based upon the pioneering paper of Wilson (1983). We use his 
formalism, but extended and with a more numerically stable normalization scheme (from Abbott & Schaefer 1986, 
who presented the first calculation of the large angle anisotropics in a closed universe). This normalization is also 
used by Hu & Sugiyama (1995). We have also drawn on the work of Gouda, Sugiyama & Sasaki (1991) which 
provides some (but not all) of the steps missing from both Wilson (1983) and Wilson's (1981) Ph.D. thesis. Our 
intention is to provide a complete derivation of the equations which are needed to evolve cosmological perturbations 
in a curved background. 

In order to keep this appendix to a manageable size, we assume that the reader is familiar with the derivation 
and implimentation of the fluid and Boltzmann equations in a flat universe. An excellent reference for this case 
is Ma & Bertschinger (1995), and our notation is based partly upon that work. Further details can be found in 
HSSW and references therein. We will only concern ourselves with scalar perturbations, although following Abbott 
& Schaefer (1986) it would be relatively straightforward to generalize to vectors or tensors. An analytic treatment 
of the anisotropics induced by tensors has been presented by Allen, Caldwell & Koranda (1995). 

A.l The Metric 

We define our metric to be 

ds 2 = -dt 2 +a 2 (t)( 7ij + h ij )dx i dx j (26) 

= a 2 (t) [-d V 2 + (jij + hij) dx'dx j ] , (27) 

where jij is a 3-metric with constant curvature K = ffg(f2tot — 1)> V = J dt/a(t) is the conformal time and hij is the 
metric perturbation in synchronous gauge (Lifshitz 1946, Landau & Lifshitz 1975, §97, Weinberg 1972, Peebles 1980, 
Peebles 1993). In our notation latin indices run over spatial components from 1, 2, 3 and greek indices over space 
and (conformal) time components 0, 1, 2, 3. 

By definition, the geodesic equations are those which extremize the Lagrangian, C = \g IJlv x^x 1 ' (where an overdot 
here only represents a derivative with respect to some parameter, A, along the path). Equating the Euler-Lagrange 
equations for C with the other form of the geodesic equation, 

rl 2 ^ „ dx v dxP , . 

r ^-7T-7r = °> ( 28 ) 



d\ 2 vp d\ d\ 
we can read off the connection coefficients: 

T° =(h/a) 

T °j = («/«) (lij + h ij) + \hij 

r>,. =(d/a)6l + \hl ( 29 ) 

ri, = (3) ^ + ^{iK)\ k + (Kh + (h^)- 

Here ( 3 )r is the connection for the spatial metric jij of constant curvature, and a subscript \i means a covariant 
derivative wth respect to this metric. This is distinguished from ; a which refers to a covariant derivative with 
respect to the full 4-dimensional metric. An overdot means d/drj henceforth. 
All that we will need explicitly from jij is the curvature tensor 

^R™ k = K(S™ 7lk -S™ 7lJ ), (30) 

which allows us to commute covariant derivatives, e.g. in the simplest case Viy^ — = ^-R™^™- However, we 
shall be concerned with commuting pairs of indices among arbitrarily many: 

V\i 1 i 2 ...i l ij...k — (V\i 1 i 2 ---iti) \j...k 

= • y, i:i>i! < <.< r < -.' •••<■ ) • ( 31 ) 

V ™=i / \j—k 
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Note that indices after the pair to be commuted go along for the ride. 
For deriving the field equations we will also need the Ricci tensor: 

Rap = (r^), 7 - (TZ 7 \p + Tl 7 T^ - Tj a T° 7 (32) 

A. 2 Basis Functions 

In linear theory, modes with different k evolve independently, so working in a wavenumber basis is advantageous. 
Although in a non-flat universe the Fourier transform is not easily defined, we can still work in terms of the 

eigenfunctions of the Laplacian operator, which we shall call Q(x). This is just e tk ' x in flat space; a detailed 
discussion of representations and spectra of Q(x) can be found in Abbott & Schaefer (1986). We will only need 
the defining relation: 

V 2 Q = t'JQi,-,- = -k 2 Q 

= -(P 2 -1)KQ. (33) 

The last equality defines /3 for fl t ot > 1- I n this case the spectrum of these functions are the integers /3 = 3, 4, 5 . . ., 
since the spatial hypersurfaces have a spherical geometry (cf. Fourier series on a circle). The modes /3 = 1,2 
correspond to pure gauge modes (Abbott & Schaefer 1986), and are not included in the spectrum. 

We can expand our perturbed quantities in terms of these functions and their covariant derivatives. For example 
the metric perturbation can be written 

hij(x, T]) = ^[h + H] Qfij + k- 2 HQ\ij, (34) 
where h is the trace of the perturbation and H is the traceless part. 

A. 3 Einstein and Fluid Equations 

The evolution equations for the background metric and its perturbation follow from Einstein's equations, which 
can be written 

R al) = 8ttG (r a p - ^g a pTj = 87rGT a/3 , (35) 

where T a f} is the stress-energy tensor and T is its trace. The background stress-energy tensor is of the perfect fluid 
form (with density p and pressure p) plus a perturbation: 

T ° = -p(i + Q6), (36) 

= -(P + P)v(ik)- 1 Q U , (37) 

Tj = (p+Q[«p])«j-+Sj, (38) 

with T,'j the anisotropic stress (E,- = 0), which we will not need explicitly. 

The required evolution equations come from the time-time and time-space Einstein's equations. The unperturbed 
equations are the familiar Friedmann-Robertson- Walker equations for the evolution of the scale factor a(r/). To 
first order in the perturbation (first-order quantities denoted by the pre-superscript (1)): 

(1) Roo = -\(^+^Q (39) 

and «T 00 = ^a 2 (pS + 3[6p])Q, (40) 

so that 

4- (ah) = 8irGa 3 p(l + 3w)6, (41) 
ar\ V / 

with w = pi p. The time-space equation involves 

2x MRU = (h)) -Ti(h 

"1 



J/ \j v >\i 



-hQ\ 
2 



- r + H) Q|,. + k- 2 H(Q\% 
-l{h + Hk- 2 (k 2 -Slir^Qi,-. (42) 
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In the last line we have made use of the following identity (coming from equation 31): 

7 jk Q\ijk = 7 jk Q\jik 



Q\jki + ^RfikQ\r 



= ~k 2 Q\, + 7 3k K (S™ 7jk - S? 7ji ) Q ]m 

= (-k 2 + 2K) Q\i. (43) 

The stress energy tensor is 

Wfoi = -a 2 (p + p)v(ik)- 1 Q li . (44) 
If we redefine H to absorb the factor (1 — 3K/k 2 ) in (42), then the evolution equation is 

k 2 /■ ■ \ 

y (h + Hj = -8irGa 2 p(l + w)6, (45) 

where 6 = ik ■ v. 

The equations for the non-relativistic fluids (CDM and baryons) can be simply obtained from the conservation 
equation for the stress-energy tensor: = 0, or 



(V=9T£) = -a.,-.,,'! -. (46) 



We will show later that for fluids without anisotropic stress the fluid equations are unchanged from their flat 
space counterparts. Explicit evaluation of these equations in the K = 0, p = case is straightforward. Using 
yj—g = a 4 (l + \h), we find for the a = equation 

s = -(e + hy (47) 

while with a = i we obtain 

= ~6. (48) 

A. 4 Legendre Tensors 

The radiation distribution function will exhibit an azimuthal symmetry, so we would like to use an expanion 
in Legendre polynomials. To avoid introducing a fixed basis, it is convenient to define the following functions of 
position x and an angle h, 

P(o) = 1, 

^i) = n\ 

P(2) = ^(3n'V- 7 <J ') 

Pg* = i (5n*Vn* - 3^ ij n k ^ , 

(2l+l)n^Pi 2 l) -- ieie+l) = ^7 (<1<a P( / ^) +l) + (^+l)P( / +'i) +1 , (49) 

where parentheses denote symmetrization in the indices. Each is a function of position through jij only (so its 
covariant derivative vanishes), and is explicitly symmetric. As can be seen from the definition, is a "Legendre 
polynomial with dangling indices", i.e. rather than a polynomial of a dot-product of vectors, the tensor has a 
free index for each power of the vector argument. Later this will allow us to perform an expansion in Legendre 
polynomials by contracting tensors in an explicitly coordinate free way. Clearly 

Pfp =Pi(k-n) (50) 

for any unit vector k. 
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Some properties of P^ which will be useful are 

JijP, 



ij---k 



]*(() 



tup; 



ij---k 



□J-"* 
■(/) ' 



which are now proven by induction. Both are true for £ = 0, 1, 2, 3 by inspection. 
Assuming (52) holds true for all P^ with i < £ + 1, we have 



= 7TT [{(2*-i)« (<1 p£:^ 
= 7^[(«+i)» (, ' 1 ^-5 ) -^p] 

where the factors — l)/2 and £ are the number of terms in the symmetric products. Also 

2£+ 1 



Putting these together gives equation (52). 

Similarly assuming (51) holds true for all P^ with i < £ + 1, we have 



(f + 1W • p 8 i--- 8 <+i 

^ -r ±) /iii 2 - r (£_|_i) 



(2£ + ih.- 1 .-y*' i pjp +i) - ^^P^p^p 



(2^+l) ?p .-3-w +1 2 

(£+1) £(£+ 1)" V^^- 1 ) 



3P 



2(£-i)p--; 



«3"-«* + l 

(£-1) 



ppp [2(2^+1) -2(2^ + 1) 



0, 



where the factors of £ + 1 and £(£ + l)/2 are the number of terms in the symmetric products. 

A. 5 Recurrence Relation 

Later we will expand our radiation distribution function in terms of mode functions 

G t = (-k)- 1 Q\ il ... it (x-j)P$" it (x,h). 

The interesting property of these functions for our purpose is the recurrence relation that they satisfy, 

+ 1 



n l G t \i = k 



2£+l 



2£+l 



Gi+i 



where 



(51) 
(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



1 - (£ 2 - l)K/k 2 £>1. 



(58) 



Notice that in terms of (3 we can write ki = \J(P 2 — £ 2 )/(j3 2 — 1) for £ > 1. 

The proof of this relation is simply an exercise in combinatorics and commuting covariant derivatives using (31). 
Let us first prove some identities which will be crucial. The simplest is 



f)i ■ ■ ■ r>' p'i—'t — f),- 



l (i P, 



if-ii) 



-^-W-kf-'G^ 



(59) 
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which comes from commuting the i to all £ + 1 possible positions (recall that P^ is symmetric in all indices). An 
index swap of i in the (n + l)th place with i„ in the nth place gives (n — 1) factors of KGi-\. The number of 
factors obtained when symmetrizing n' P(t) is thus 

Ei(i + 1) = ^ 2 -1), (60) 
i=o 

which when divided by the £ + 1 terms gives the relation above. 
A more taxing but just as conceptually simple identity is 

n >ip'>-M 2 ( , 2 1(1+1) , (1 + 2)! \ t _ x 

which comes from looking at contractions of the two indices, i and ii, in arbitrary positions on Q| with y tte P^_^ _1 . 

If the indices are in the first two positions then one obtains (— k 2 Q)\i 1 ...i e _ 1 from the eigenvalue equation for Q. 
Thus the problem is: how many factors of K does commuting indices introduce before these two indices are in the 
first two positions? If we decide that ii will always come before i (otherwise use the symmetry of j"' and relabel 
these indices) then we first shift ii from the nth position to the first, then shift i from the mth position to the 
second position. Moving something from the nth position to the (n — l)th position, assuming that it is the first 
index on j"' gives (n — 2) factors of K. Moving this index all the way to the left thus introduces (n — l)(n — 2)/2 
factors of K. The second index is more complicated since now the possible contractions of indices allow jijj' J = 3 
and 7ij7 jfc = Sf . One finds that moving the mth index to the second position gives m(m — l)/2 — 1 factors of K. 
The total number of K , s is thus 



^ (n-l)(n-2) 
^ 2 

1 m=n+l 



m(m — 1) 



V ' (62) 



6(^-2)! 



Now we are in a position to derive the recurrence relation: 



-*)-' 








-*)-' 




> + ±£(£-l)K(-k) l - 1 G i - 1 




-*)-' 


l -£(£ - i)K(-k f 


l ( 1 ( 

Gt-i + Qii.-.iti y 2 £+ i 7 




k 1+1 G l+1 + ( k)~ l 


lt(£-l)K(-k) t - 1 G l - 1 + 





\ £ + 1 
+ 2£+l 



2£ +!£(£+ 1) 



2 K J 6(£-2)l ) K ' 



2£+l 



K/Gf 



2£+l 



(63) 



Let us try to understand this relation physically. As we shall see in the next section this recurrence relation 
describes the way power is transferred between £ modes by propagation. Consider first the case r^tot — 1? or 
K = 0. A given physical scale subtends a smaller and smaller angle (higher £) as the distance between emission and 
reception of the photons increases. This is characterized by the coupling of the £ modes in the recurrence relation, 
which allow power to be transferred from low £ to high £. In the case of a non-flat universe (K ^ 0), the manner 
in which power is transferred from low to high £ is modified by the effects of curvature (geodesies are no longer 
straight lines). For an open universe (K < 0) the coupling between modes is enhanced, making it easier to transfer 
power to higher £ (smaller scales). This is because geodesies in an open universe diverge, thus a fixed physical 
scale subtends a smaller angular scale at the same distance in an open universe. In a closed universe (K > 0) the 
coupling between modes is decreased, making it harder to move power to higher £. In fact when k 2 / K = £ 2 — 1, 
m = and there is no coupling to higher modes! In a closed universe there is a finite heirarchy for each mode. This 
reflects the fact that in a closed universe, there is a maximum physical scale that is probed by geodesies with fixed 
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angular separation. Note also that once power reaches this I it will reflect back to lower £, showing that the same 
physical scale begins to subtend a larger angle after the maximum separation of the geodesies has been reached 
(at the "equator"). It is this property which is reponsible for the behaviour of the position of the first peak in the 
power spectrum (see Fig. 3). 

One last important property of the Gi is their normalization. Notice that Q(x; k) depends on the direction of 
k. Recalling that is a symmetric tensor and (50, 51), we have that 



dSl k G t oc J dSl k QPi(k ■ h). (64) 
To fix the normalization of these functions we use the recurrence relation and 

dn k I ^d 3 x [G/G/_ = 0, (65) 



to find 



A. 6 The Boltzmann Equation 



j dSl k \\G\\\ = K f j dSl k IIGJUII. (66) 



We are now in a position to derive the Boltzmann equation for the radiation distribution function, which is 
a function of the 6 dimensional phase space for the radiation (x, p) and the conformal time rj. We scale out the 
redshifting of the photon momentum, p, by the definition q = ap and write q = qri. Integrating f(x, h, q, rj) over q 
gives the brightness perturbation, F(x, h, rj), which we can expand in terms of our basis functions as 



F(x,h, V ) = J2J2 F ^^nIl Kl ) G ^,n), (67) 
k I \ I ) 

where the factors of ki have been included to provide a convenient normalization of the functions (see equation 66). 
By including the factors of ki we ensure that the expression for Ci in terms of Ft will be formally unchanged from 
the K = case. Absorbing the very large factorials involved in Yii K t ' ls a l so important to ensure numerical stability 
in the computation of the radiation power spectrum. 

The Boltzmann equation for the radiation can be written schematically as 

DF _ dF d .dF 

Drj drj dr\ ^ ' ^ dq 

dF -dF . t dF .dF 
~ -t— + Xt-^ + n*—^ + q — 
or] ox* on* oq 

= Collision Term. (68) 

For the photons the RHS is the collision integral for Thomson scattering, which can be found in Hu, Scott & Silk 
(1994) and Dodelson & Jubas (1995), while for the neutrinos the RHS is zero. In the second line we have expanded 
the rj derivative, which is formally the change along the photon path, n l F\i, and written it schematically in terms 
of x' and h' , hence the approximately equal sign. In a flat universe the h' term vanishes: the geodesic is a straight 
line. It is this term which gives rise the the ki factors in the recurrence relation for the Gi. 
We can evalulate the q term from the geodesic equation: 

? ; i j 

q = -- hijn n J 

= ( t(-hQ + 2H{l-ZK/k 2 )- 1 G 2 y (69) 

Putting these terms together, and using the recurrence relation (57) to rewrite the derivative along the path, we 
find the following (recall 6 = ik ■ v): 

4 2 • 

6 = — 6--h 
3 3 
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k 2 (-S - ^ 2 F 2 ) + an e a T (0 B - 0) 



F 2 
Ft 



3 4 

k 2 - -kn 3 F 3 + -k 2 1 H 
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an e <JTF 2 



21- 



21- 



-Kl + lFl + l 



an e a T Ft, 



(70) 



where for the neutrinos the terms proportional to cry, coming from the collision term, are absent. The 1 = 2 mode 
is related to the anisotropic stress. As can be seen, for I < 1 the equations are formally identical to those for K = 0. 

The modification of the equations for photon polarization (Bond & Efstathiou 1984, Kosowsky 1995, Zaldarriaga 
& Harari 1995) simply requires inserting factors of ki in the mode coupling terms just as above. Since for S7 to t > 1 
the sound horizon can subtend a larger angular scale than in a flat universe, the polarization-to-temperature ratio 
can be larger in a closed universe, as one would expect by analogy with the open case (Zaldarriaga & Harari 1995). 

The remaining (perturbed Einstein and fluid) equations are gathered here for completeness: 



drj 



ah 



&CDM 

Sb 

OfJDM 

B 



(87rG)a 3 ^ft(l + 3 Wi )^ 

i 

-{^G)a 2 Y,Pi( l + w i) e i 



-h 

2 

-0B ~ \h 



--Ob + an eVT (0 7 - Ob)- 
a 3p B 



(71) 



where the sum J^. is over the components: CDM, baryons, photons, neutrinos. 

The initial conditions for the metric, density and velocity perturbations are unchanged by curvature, as one would 
expect physically, and can be found in Ma & Bertschinger (1995). The initial value of the neutrino anisotropic 
stress acquires a factor of k 2 in a closed universe. Due to our choice of normalization the expression for Ci is the 
same as for the flat space case: 



(2£ + l) 2 C t 



VK 
~8T 



oo 



13=3 



,(/3) ^\F t (p, m )Y 



(72) 



where r/o is the conformal time today. The factor of yK multiplying the sum over (3 makes the limit yK — > 
smooth, and should also appear in the definition of the moments of the matter power spectrum. 

The primordial power spectrum, P pr ; m (/3), comes from a theory of fluctuation generation. Rather than perform 
calculations of the primordial spectrum for specific inflationary potentials, we have assumed that the spectrum 
is either a power-law in /3 (with spectral index unity), or that it corresponds to constant fluctuations in the 
gravitational potential $ per logarithmic interval in wavenumber (which in the flat case is the same as the previous 
assumption). For this latter case, generalizing (Lyth & Stewart 1990, or see White & Bunn 1995) gives 

(B 2 - 4) 2 

which we shall refer to as "scale invariant". 
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